NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



OR /<&s4$ 


( 


A CORRELATION METHOD TO PREDICT THE SURFACE 
FRESSURE DISTRIBUTION ON AN INFINITE PLATE 
OR A BODY OF REVOLUTION FROM WHICH 
A JET IS ISSUING 


by 


Stanley c. Perkins, Jr. and 
Michael R. Mendenhall 


(NASA-CR-152345) A CORRELATION METHOD TO N80-32339 

PREDICT THE SURFACE PRESSURE DISTRIBUTION OF 
AN INFINITE PLATE OR A BODY OF REVOLUTION 

FROM WHICH A JET IS ISSUING Final Report, 1 Uiiclas 

Dec. 1978 (Nielsen Engineering and Research, G3/02 34392 


NEAR TR 211 

\ 

1 ' 

January 1980 / 


Prepared under Contract No. NAS2-10125 

for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Ames Research Center 

by 

NIELSEN ENGINEERING & RESEARCH, INC. 

510 Clyde Avenue, Mountain View, CA 94043 
Telephone (415) 968-9457 



1, Report No, 2. Government Accession No, 

NASA CR-152,345 

3, Recipient’s Catalog No, 

4. Title and Subtitle A CORRELATION METHOD TO PREDICT 
THE SURFACE PRESSURE DISTRIBUTION ON AN 
INFINITE PLATE OR A BODY OF REVOLUTION FROM 
WHICH A JET IS ISSUING 

5. Report Date 

May 1980 

6. Performing Organization Code 

549/C 

7, Author(i) 

Stanley C. Perkins, Jr. and Michael R. 
Mendenhall 

8. Performing Organization Report No, 

NEAR TR 211 

10, Wch»j Unit No, 
\ 

9, Performing Organization Name and Address 

Nielsen Engineering & Research, Inc. 
510 Clyde Avenue 
Mountain View, CA 94043 

11. Contract or Grant No. 

NAS 2-10125 

13, Type p ® nt * *’ er *°d Covered 

12/7/78 - 5/1/80 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administratio; 
Ames Research Center 
Moffett Field, CA 94035 

1 14, Sponsoring Agency Code 

15. Supplementary Notes 

Ames Technical Monitor, K. Aoyagi 


is. Abstract A correlation method to predict pressures induced on an 
infinite plate by a jet exhausting normal to the plate into a 
subsonic free stream has been extended to jets exhausting at 
angles to the plate and to jets exhausting normal to the surface 
of a body of revolution. The complete method consists of an 
analytical method which models the blockage and entrainment 
properties of the jet and an empirical correlation which accounts 
for viscous effects. For the flat plate case, the method is 
applicable to jet velocity ratios up to ten, jet inclination 
angles up to 45° from the normal, and radial distances up to five 
diameters from the jet. For the body of revolution case, the 
method is applicable to a body at zero degrees angle of attack, 
jet velocity ratios 1.96 and 3.43, circumferential angles around 
the body up to 25° from the jet, axial distances up to seven 
diameters from the jet, and jet-to-body diameter ratios less than 
0.1. Correlation curves and data comparisons are presented for 
inclined jets issuing from a flat plate with velocity ratios one 
to twelve. Correlation curves are presented for jets issuing 
normally from the surface of a body of revolution for velocity 
ratios 1.96 and 3.43. 


17. Key Words (Suggested by Author(si) 

Subsonic Crossflow 
Jet in a Crossflow 
Pressure Distribution 
V/STOL 

18. Distribution Statement 

Unclassified - Unlimited 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

198 

22, Price* 


‘For sale by the National Technical Information Service, Springfield, Virginia 221 61 





TABLE OP CONTENTS 


Section 
SUMMARY 
INTRODUCTION 
LIST OF SYMBOLS 
APPROACH 
Flat Plate 
Jet model 
Correlation 
Body of Revolution 
Body model 
Jet model 
Correlation 
RESULTS 

Flat Plate 
Body of Revolution 
CONCLUSIONS 
RECOMMENDATIONS 
REFERENCES 

FIGURES 1 THROUGH 36 


ill 


Jl 




A CORRELATION METHOD TO PREDICT THE SURFACE 
PRESSURE DISTRIBUTION ON AN INFINITE PLATE 
OR A BODY OF REVOLUTION FROM WHICH 
A JET IS ISSUING 


by 


Stanley C. Perkins, Jr. and Michael R. Mendenhall 


SUMMARY 

A correlation metho.d to predict pressures induced on an 
infinite plate by a jet exhausting normal to the plate into a 
subsonic free stream has been extended to jets exhausting at 
angles to the plate and to jets exhausting normal to the surface 
of a body of revolution. The complete method consists of an 
analytical method which models the blockage and entrainment 
properties of the jet and an empirical correlation which 
accounts for viscous effects. For the flat plate case, the 
method is applicable to jet velocity ratios up to ten, jet 
inclination angles up to 45° from the normal, and radial 
distances up to five diameters from the jet. For the body of 
revolution case, the method is applicable to a body at zero 
degrees angle of attack, jet velocity ratios 1.96 and 3.43, 
circumferential angles around the body up to 25° from the 
jet, axial distances up to seven diameters from the jet, and 
jet-to-body diameter ratios less than 0.1. Correlation curves 
and data comparisons are presented for inclined jets issuing 
from a flat plate with velocity ratios one to twelve. 

Correlation curves are presented for jets issuing normally from 
the surface of a body of revolution for velocity ratios 1.96 and 
3.43. 






INTRODUCTION 


Increased interest in V/STOL configurations which utilize 
lift-jet engines mounted in the wing, pod, or the fuselage has 
prompted experimental and analytical studies of jet 
interference effects on infinite flat plates and bodies of 
revolution. Experimental studies on infinite flat plates (refs. 
1-7) and bodies of revolution (refs. 3, and 8-12) have shown 
that the jet produces a region of positive pressure upstream of 
the jet and regions of negative pressure downstream and to the 
side of the jet. These studies also show a very strong viscous 
region immediately downstream of the jet. It is in this 
viscous-dominated region that inviscid analytic methods, such as 
those presented in references 13-17, have not been successful in 
predicting surface pressures. As an alternative to purely 
inviscid methods, an empirical method to predict surface 
pressures on a flat plate in a crossflow with a jet exhausting 
normal to the surface was developed by Nielsen Engineering & 
Research, Inc. (NEAR) and is presented in reference 18. This 
method consists of an analytical jet model for the inviscid part 
and empirically-derived factors to account for viscous effects. 
These factors, which are obtained from a correlation of the 
difference between analytically predicted surface pressures and 
measured results, are presented as a function of jet velocity 
ratio and position on the plate . 

In this report, results of a two-phase study to extend the 
empirical method of reference 18 are presented. The first phase 
deals with jets exhausting at an angle to an infinite plate. 

The analytical jet model of reference 18 is used in conjunction 
with available experimental data to obtain empirical correlation 
factors. These viscous correlation factors are presented for a 
range of jet inclination angles and jet velocity ratios, and 


predicted pressure distributions are compared with experi- 
mental data from several sources. Also included as a part of 
the first phase are comparisons of new normal- jet data from 
NASA/Ames Research Center and results obtained using the 
prediction method with the viscous correlation factors. 

The second phase deals with jets exhausting normal to the 
surface of a body of revolution. The analytical model utilizes 
the NSRDC version of the Douglas-Neumann program (refs. 19 and 
20) to model the body/ jet-surface mutual interference effects 
and the sink distribution model of reference 14 to represent 
entrainment effects. The jet boundary and position are 
determined from empirical observations. Development of viscous 
correlation factors for jets exhausting from a body of 
revolution are presented for a limited range of jet velocity 
ratios . Recommendations for improving the present correlation 
method for jets exhausting from a body of revolution are 
discussed . 








LIST OP SYMBOLS 


pressure coefficient, p - p^/q^ 
jet diameter at the exit plane 
unit normal at- a point, see figure 4 
static pressure 

1 2 

dynamic pressure, •jPcqV 

radial distance along the plate from the center 
of the jet to any field point on the plate 

jet radius at the exit plane 

curve length of the jet axis 

curve length of potential core 

local jet radius 

velocity components in the x,y, and z directions, 
respectively 

jet Velocity at the exit plane 

constant free-stream velocity 

plate coordinate system fixed at the center of 
the jet exit plane, positive x is upstream 

jet coordinate system fixed at center of the jet 
exit plane, positive x^ is downstream 

polar angle, measured clockwise from the positive 
x-axis in the plate x-y plane, see figure 1 

correlation increment of pressure coefficient 

initial inclination angle of jet centerline, 
measured from the positive x-axis in the x-z 
plane, <Sj = 90°-8; see figure 1 

initial inclination angle of jet centerline, 
measured from the positive zj axis in the Xj-Zj 
plane, 0=0° for a jet issuing normal to tne 
free stream, see figure 1 
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LIST OF SYMBOLS (Concluded) 


p OT density 

<f> circumferential angle/ measured from top of body 

of revolution, see figure 17 


Subscripts 

j jet quantity 

max maximum value 

oo free-stream quantity 
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APPROACH 


The problem of interest is the prediction of the surface 
pressure distribution in the vicinity of a jet issuing from a 
flat or curved surface in the presence of a free-stream flow 
along the surface. The flat plate case is illustrated in figure 
1. If one uses a potential flow model to describe the flow, 
including jet entrainment, blockage and curvature, one can 
predict the surface pressures upstream of the jet (3 less than 
60°) reasonably well, but the discrepancies become 
increasingly large as 3 increases beyond 90°. The source of 
the differences is a viscous interaction between the flow around 
the jet and the surface flow, which, cannot be modeled properly 
with present analytical methods due to lack of detailed flow 
data in the interaction region. 


In the work of reference 18, an approach was adopted for the 
flat plate case in which a potential model was selected which 
gave good agreement with surface pressures upstream of the jet, 
and the differences between the predicted and measured pressures 
around the jet were correlated into a viscous correction. 
Accordingly, the pressure coefficient at every point on the 
plate can be expressed as 


C = C 
P P 


+ AG 

potential P 


viscous 


( 1 ) 


where ^Cpjviscous the correlation factor which represents 
viscous or non-potential effects. The potential portion of the 
above equation is calculated using the Bernoulli equation in the 
form 
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( 2 ) 


( (V + u) 2 + v 2 ) 

°’ = 1 -f * ( 

In the current work# the same approach is taken. The first 
task, then, is to establish a potential jet model that agrees 
well with the upstream pressures, considering first the case of 
the flat plate with the jet issuing at an angle <$ • other than 
90° and secondly the case of a jet issuing normally from the 
surface of a body of revolution. After the model is selected 
and verified, the viscous correlations are obtained. In the 
following sections, the development of the jet model and 
correlation for the two cases is described. 


Plat Plate 

A sketch of an expanding jet of initial velocity Vj and 
initial inclination angle 0 emerging from an infinite plate into 
a subsonic crossflow of velocity is shown in figure 1. The 
flat plate method presented herein is applicable to low speed 
crossflow velocities and jet Mach numbers less than one. 

Jet model .- Details concerning the specification of the jet 
centerline path, the blockage and entrainment models, and the 
spreading rate of a jet exhausting normal to a flat plate are 
presented in reference 18. The same basic jet model, with some 
modifications made to the blockage and entrainment models, is 
used for a jet exhausting at an angle to a flat plate. A 
general description of the various components of the jet model 
is presented in the following paragraphs. 


7 


The path of the jet centerline is specified using the 
empirical relation developed by Margason (ref. 21); 
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Figure 2 shows jet centerline shapes given by equation (3) for a 
range of velocity ratios with <$j = 75°, 60°, 45°, and 30°. 

As shown in this figure, the spread in centerline position 
with velocity ratio is reduced as the jet becomes more parallel 
to the crossflow. 


A comparison of measured and predicted jet centerline shapes 

is shown in figure 3 for 7 . /v = 12 and 0 = 0°, 15°, 30°, 45°, 

_ J 

and 60°. The good agreement between the measured centerline 
shapes, which were determined by locating the position of 
maximum total pressure in successive yj - z.. planes behind the 
jet, and those obtained using equation (3) lends confidence to 
the use of this equation for obtaining centerline shapes fox' 
jets in a cro/ssflow. 

Blockage effects of the jet are modeled by representing the 
specified wake surface with a series of vortex quadrilateral 
panels on which a boundary condition of flow tangency due to the 
free stream is enforced. The vortex quadrilateral panel 
blockage model described herein is based on a method developed 
by Maskew (ref. 22). 

The boundary of the jet wake, specified by the centerline 
path and the wake radius along the centerline, is divided into 
length segments as illustrated in figure 4. Each length segment 
is further divided into a number of circumferential segments. A 
slice through the wake, normal to the centerline, shows the 
circular cross section wake to be represented by a series of 


8 


straight segments. The control point on each panel is located 
t.fc -the centroid of the panel, and the vector n is the unit 
normal to the panel, positive when directed outward. The 
boundary condition satisfied at each control point is 

V • n 0 (4) 


Where V is the vector sum of the free-stream velocity and the 
velocity induced at a panel control point by all vortex 
quadrilateral panels making up the blockage model. 

Implementation of this boundary condition results in a set of 
linear simultaneous equations which can be solved for the 

unknown vortex quadrilateral strengths. 

* 

After obtaining the vortex strengths for each panel, the 
field point velocity induced by a single vortex quadrilateral is 
obtained by summing the contribution from each of the four sides 
of the quadrilateral. A summation over all the panels 
representing the jet surface produces the velocity induced at a 
point by the complete blockage model. 

As a jet is inclined from the normal to the plate, the shape 
of the jet exit at the plate becomes elliptical. In order to 
model this shape, modification of the panel layout in the region 
near the plate is necessary. Figure 5 shows a side view of the 
panel layout of a jet in the region near the plate for Vj/v^ = 
12.0, 0 = 45°. The "transition" disc between the elliptical 
disc at the plate and the first circular disc is required 
to maintain panels of reasonable aspect ratio (height/width) in 
the front part of the jet. A disc, as indicated in the figure, 
is a slice through the jet cross-section at the end of a row of 
panels. Each disc following the third disc is circular. 
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It is roted that the predicted plate pressures in the region 
immediately behind the jet become very sensitive to panel 
length for large jet inclination angles. This is a numerical 
problem due to end effects from vortex quadrilateral panels 
which come closer to the plate in the region behind the jet as 
the inclination angle increases. It is necessary that low 
aspect ratio panels be used for jets exhausting at large angles 
to a flat plate. 

The spreading rate for a jet exhausting from a flat plate is 
calculated using a method described in reference 18. This 
method utilizes the empirical method of reference 23 and data 
from reference 24. The empirical method of Bowley and Sucec 
(ref. 23) is directly applicable to a jet exhausting at an 
angle. It is assumed, due to lack of data, that the loc-^1 
centerline velocity data of reference 24 for jets exhau- Lng 
normal to the plate can be used for jets exhausting at an angle. 
Total pressure profiles presented in reference 7 for a jet 
exhausting at a range of inclination angles to a plate indicate 
that this assumption is valid. 

Assumptions are also made, as discussed in detail in 
reference 18, that the jet is axisymmetric and the potential 
core length is three jet diameters. It is generally accepted 
that potential core length is a function of jet velocity ratio 
and inclination angle; however, the variation of 
experimentally-determined potential core length with these 
parameters, a;; shown in figure 6, is considerable for different 
tests and there is insufficient data on which to base a 
functional relation. Consequently, a potential core length of 
three jet diameters was selected as a reasonable average length 
based on the range of data shown. 
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Figure 7 presents a series of expansion rate curves which 
show nondimensional jet radius (t/D) versus nondiraensional 
distance along the jet center]. ine (S a /R) . These curves are used 
to represent the expansion of the circular cross section of the 
jet (ref. 18). Some comments regarding the curves in figure 7 
follow. First, data from reference 24 were used to obtain 
expansion curves for Vj/V M =4.0, 6.0, and 8.0. All other 
expansion curves were obtained by interpolation or extrapolation 
of these curves. Second, the Vj/Vco = 12.0 expansion rate 
curve is not shown, since the expansion curves for Vj/Voo = 

10.0 and 12.0 are essentially identical. 

Figure 8 presents plots of (t/D) max versus V j/V m , where 
(t/D) max is the value o'f t/D at S a /R = 20. This value of jet 
centerline length (S a /R) is used as the upper integration limit 
for the jet entrainment model and represents the total length of 
the jet modeled. It was found (ref. 18) that increasing the 
length of the jet beyond 20 initial jet radii increased 
computation time without affecting change in the final results. 
Choosing S a /R = 20 as the length of the jet centerline and 
assuming (t/D) max + 1*0 as Vj/V^ + 0, a curve is faired 
through this point and the values of (t/D) for V-i/V- = 

4.0, 6.0, and 8.0. These plots are then used when interpolating 
or extrapolating the Vj/V^ =4.0, 6.0, and 8.0 expansion 
curves to obtain a spreading rate for arbitrary Vj/V^ (see 
ref . 18 ) . 

Jet entrainment effects for a jet exhausting at an angle to 
a plate are obtained using a modified version of an empirical 
method developed by Yeh (ref. 14). Experimental data from 
reference 7 indicate that the entrainment rate of a jet is 
lowered as jet inclination angle (0) is increased; however, the 
amount by which entrainment is reduced is not known. Since 
Yeh's method was developed for jets exhausting normal to a flat 
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plate, it was decided to use an effective velocity ratio given 
by the normal component of the initial jet velocity ratio. This 
relationship is given by: 




effective 


V i 

cos 0 

V oo 


(5) 


where 0 is the initial inclination angle of the jet measured 
from the normal to the plate (see fig. 1). 

Predicted results were obtained using these entrainment and 
blockage models, and comparisons with experimental data over a 
wide range of velocity ratios and inclination angles were made. 
These data comparisons, shown in figures 9 through 12, are for 
the region of the plate where viscous effects are believed to be 
small (0° i p < 60°) based on experimental observations 
(ref. 5). Data comparisons in this region are a true test of 
the ability of the potential flow model to accurately predict 
pressures on the plate in the absence of large viscous effects. 

AS can be seen in these figures, varying degrees of success 
were obtained using the reduced entrainment model. In general, 
predicted results compare well with experimental results for low 
velocity ratios (V^/V < 7) and poorly with experimental 

results for high velocity ratios, especially for 0 =* 15° and 
45°. For cases in, which good agreement between theory and 
experiment were obtained, differences in the trends exhibited by 
the data comparisons make it difficult to assess the reduced 
entrainment model. For example, data comparisons for Vj/V^, = 
4.0, 0 = 30° [fig. 10(b)], indicate good agreement between 
the data of reference 6 and predicted results obtained using the 
standard entrainment model. However, this same figure shows 
that further reduction of the effective velocity ratio is needed 
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to improve agreement between predicted results and data from 
reference 28 for 8=0° and 30°. Differences between 
data sets and between data and theory, as seen in this figure, 
may be due to real jet effects, such as initial velocity 
profile, turbulence level, and potential core length which are 
not included in the present jet model. 

Other low velocity ratio cases were found in which the 
differences between predicted and experimental results may be 
due to real jet effects or plate edge effects. For example, 
figure 13 shows a comparison between measured and predicted 
pressures for Vj/V^ = 6.67, 8 = 15°, 8 = 60°. The increment 
between theory and experiment at r/D = 5.0 is AC^ = .074. 
Applying this increment to the theory improves the comparison 
between theory and experiment over the entire r/D range. This 
phenomena may be due to edge effects or real jet effects, since 
the increment is usually constant in the region of the plate 
where viscous effects are small (8 = 60°). Since such 
effects probably vary with distance from the jet or edge of the 
plate, the increment probably varies with r/D and cannot simply 
be added to every theoretical value, Further investigation of 
this phenomena is necessary before applying such techniques to 
the theory. 
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In addition to poor agreement between theory and experiment 
for high velocity ratios, it was found that in many cases the 
differences between data sets exceeded the differences between 
theory and experiment. For example, figure 9(f) presents a 
comparison of predicted and experimental results from 
references 6 and'28 for V. /V • = 10.0, 0 = 15°. As this 

j 9 00 

figure shows, the predicted results split the difference between 
the two sets of data. 
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The disagreement between different sets of data at higher 
velocity ratios suggested the need for further examination of 
the available data. To do this, data from several sources were 
plotted in a systematic fashion to allow direct comparisons of 
trends between data sets. Figures 14 and 15 show the effects of 
inclination angle at nearly constar t velocity ratio and the 
effects of velocity ratio at constant inclination angle, 
respectively, on experimental pressure distributions at 3 = 0 °. 
Examination of these figures and of other similar data 
plots produced the results discussed in the following 
paragraphs . 

First, as exemplified in figures 14(b) and (c), the high 
velocity ratio data of reference 28 indicate more effect of 
inclination angle than is shown by the data from references 6 
and 7. Second, the high velocity ratio data of reference 28 for 
0 = 15° and 45° are nearly identical in the region upstream 
of the jet. This second result illustrates the reason 
for the large differences between data from references 6 and 28 
at 0 = 15°. Note that the data of reference 28 are still 
undergoing examination and should be considered preliminary at 
this time. 

Third, the trends exhibited by the data of references 4 and 
6 in figure 15 (i.e., decrease in pressure level with increase 
in jet velocity ratio) , are not exhibited by most of the high 
velocity ratio data of reference 28 for 0 > 15°. These same 
data from reference 28 appear to be asymptotic to a value of Cp 
other than 0 in the region upstream of the jet. Adjusting the 
pressure level of these data such that it is consistent with the 
trends exhibited by the data of references 4, 6 , and 7 (i.e., 
making it more negative and asymptotic to C =0) would greatly 
improve comparisons with predicted pressures. 
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Since data comparisons using the reduced entrainment model 

are usually poor only for cases in which there appear to be 

inconsistencies in the data, this model is used to obtain 

viscous correlation factors for jets exhausting at an angle to a 

flat plate. Based on the previous discussion of the reference 

28 data, the correlation factors for V./V >.7, 0 =15° and 45° 

j 03 

are to be considered preliminary until the apparent incon- 
sistencies in these data are resolved. 


Correlation .- The purpose of the correlation is to isolate 
the viscous or non-potential effects of the jet on the plate. 
It is expected that the predicted potential pressure dis- 
tribution away from the jet and within the range 0° £ 0 <_ 60° 
will be in good agreement with experiment. The predicted 
pressures near the jet and within the range 60° <_ 0 <_ 180° 
will likely be in poor agreement with experiment due 
to viscous effects. Assuming that the measured pressure 
distribution can be represented by a potential part and a 
viscous part, the viscous part can be determined by a 
differencing technique. 



= C 

- c 

viscous p 

experiment v 


potential 


( 6 ) 


Correlating this quantity as a function of jet velocity ratio 
and position on the plate, the predicted pressure induced by a 
jet exhausting from the plate into a crossflow is given by 



+ AC. 


potential 


viscous 


(7) 


The correlation curves presented for a jet exhausting at an 
angle to a plate were obtained in the same manner as those 
presented in reference 18. In cases where more than one data 
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set is available for a given jet velocity ratio/ jet inclination 


angle combination, (AC p ) viscous is found b Y averaging the 
correlation values obtained from each data set. In this report 
the inclination angle has been included, along with the jet 
velocity ratio, as a dominant parameter in determining pressures 
on the plate. This results in n data base which consists of a 
(ACp) array as a function of 0 and r/D for each jet velocity 
ratio and jet inclination angle. Linear interpolation is used 
to determine correlation values at any Vj/V^, 0, 0, and r/D 
values for which correlation curves have not been determined. 
These results are described in the RESULTS section. 


Body of Revolution 

The general problem of interest for the body of revolution 
case is a subsonic jet exhausting normally from the surface of a 
body at zero degrees angle of attack into an incompressible 
crossflow. The parameters of interest are jet velocity ratio, 
body shape, and jet*-to-body diameter ratio. The general method 
for the body of revolution case is applicable over a wide range 
of these parameters; however, the present method has been 
verified for only a limited range of parameters because of the 
small quantity of available data. The present method has been 
applied to a cone-cylinder model with a small jet-to-body 
diameter ratio for two low jet velocity ratios. Further 
discussion of extending the method to larger jet-to-body 
diameters and higher velocity ratios is found in following 
sections . 

Body model .- The NSRDC version of the Douglas-Neumann 
program (refs. 19 and 20) is used to model the specified body of 
revolution geometry. This program approximates the surface of 
the body by a set of plane quadrilaterals. The solution is 
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constructed in terms of a source density, which is assumed to be 
constant in each quadrilateral element. The solution satisfies 
Laplace's equation in the free-stream and a zero normal-flow 
condition on the. surface of the body. The program presently 
utilizes a maximum of 650 panels to model a three-dimensional 
body. 

As a check of the accuracy of the NSRDC program, predicted 
pressures on an ogive-cylinder body with a fineness ratio 3 nose 
at a = 10° were compared with measured data (ref. 29). The 
agreement between theory and experiment is good over the entire 
length of the body, as shown in figure 16. The difference 
between experiment and theory near <j> = 60° for x/D = 4.5 and 
6.0 may be due to viscoUs effects, which are not accounted for 
in the NSRDC program. These comparisons lend confidence to the 
use of the NSRDC method for developing viscous correlation 
factors for the body of revolution case. 

A cone-cylinder model with a jet exhausting normal to the 
surface (ref. 3) is shown in figure 17. The model was tested at 
zero angle of attack and low free-stream velocity over a range 
of jet velocity ratios. Pressures were measured on the cylinder 
surface in the vicinity of the jet. The row of pressure taps 
shown in figure 17 was rotated to obtain pressure data from the 
top of the body ( tp = 0°) to <j> = 25°. 

Due to the close proximity of the jet to the aft end of the 
body, a parametric study of the effects of modeling the 
downstream end of the body for the jet-off case was carried out. 
One model incorporated an ogive extension or boattail to close 
the body at the downstream end, and a second model incorporated 
an open end on the cylinder. Since the first model is a closed 
body, the sum of the source strengths is zero, resulting in 
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zero net flow inside the body. The second model essentially 
models the wake behind the cylinder since a net flow is produced 
due to the open end . 

Jet-off pressure distributions (fig. 18) showed a large 
effect of the ogive, such that the pressure coefficient did not 
recover to zero downstream of the cone-cylinder junction. 

Without the ogive extension, the predicted pressure coefficient 
approaches zero toward the end of the cylinder, as would be 
expected. Data comparisons for the jet-off case in figure 19 
show reasonably good agreement between the measured and 
predicted results for the model without an ogive extension. 

These comparisons indicate that the measured pressures recover 
more quickly than the predicted pressures, probably due to the 
separation at the cone-cylinder junction, which is located at 
x/D = 12.9. Such effects are not accounted for in the NSRDC 
program. As a result of this investigation, the aft end of 
the body is modeled as an open-ended cylinder. 

Jet model .- Details concerning the specification of the jet 
centerline path, the blockage and entrainment models, and the 
spreading rate of a jet exhausting normal to the surface of a 
body of revolution are presented in the following paragraphs . 

The path of the jet centerline is specified using an 
empirical relation developed by Viehweger (ref. 12). This 
relation, obtained by correlating results for a low pressure jet 
with a short potential core exhausting from a body of 
revolution, is written as 

x./D = 1.727(z./D • /q^/qT ) 2,65 (8) 
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This centerline result is based on data in the range 

16 < q./q < 290. For practical purposes, the assumption is 

made that /q“7q7 is approximately equal to V /V . , so the 

00 j Co J 

range of velocity ratios represented is 4.0 to 17.0. The jet 
centerline data from reference 12 is the only available data 
for a jet exhausting from a body of revolution. 

Figure 20 shows centerline shapes for velocity ratios 1.96 
and 3.43 as given by equation (8). Although the velocity ratios 
of interest are outside the range used to obtain equation (8), 
this equation has the advantage of having its origin in 
body-of-revolution data. Since references 3 and 12 both 
consider similar low pressure jets with short potential cores, 
it is assumed that the above centerline relation will provide 
accurate centerline shapes in the velocity ratio range of 
interest . 

Blockage effects of the jet were obtained by modeling the 
jet surface as an integral part of the cone-cylinder model. In 
this manner, the mutual interference effects between the jet 
surface and the body are accounted for simultaneously. With 
respect to the panel layout in the region of the body/ jet 
intersection, the edges of the body panels are matched to the 
edges of the jet surface panels, thereby providing a smooth 
transition from the cylinder surface to the jet surface. 

The cone-cylinder configuration (fig. 17) was modeled using 
the symmetry option of the NSRDC program, which necessitated 
modeling only one-half of the jet/body combination. A total of 
632 panels were used to model the configuration; 532 panels on 
the body and 100 panels on the jet surface. The jet is divided 
into 20 circumferential segments (10 segments per half plane) 
with panels near the body having an aspect ratio (height/width) 
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of 3.18. A fine grid on the body was utilized in the region of 
the jet in order to get better resolution of the pressures in 
that region. Predicted results could not be obtained at exactly 
<j> = 0, since results cannot be obtained along a line of symmetry 
when the symmetry option is exercised. 

Jet entrainment effects are modeled using the empirical 
method of Yeh, which is described in detail in reference 14. 
Although this method was developed for a jet exhausting from a 
flat plate, lack of an available entrainment model for a jet 
exhausting from a body of revolution necessitated the use of a 
flat plate model. The entrainment- induced velocities are 
included with the jet blockage model induced velocities in the 
boundary condition on the body surface. The entrainment induced 
velocities are not included in the boundary conditions on the 
jet surface. 

The spreading rate for a jet exhausting from a body of 
revolution is calculated using the same method used for the flat 
plate case (ref. 18). Once again, lack of spreading rate data 
or theory for a jet exhausting from a body of revolution 
necessitates the use of a flat plate method. 

It is noted here that the flat plate models for entrainment 
rate and spreading rate should give accurate results for the 
small jet-to-body diameter ratio (.064) of the cone-cylinder 
model on which data is available. The flat plate model may 
become less accurate as the jet diameter increases to an 
appreciable fraction of the body diameter, and it will be 
necessary to have additional data to evaluate differences 
between flat plate and body jets as the relative jet size 
increases . 
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Static jet dynamic-pressure decay data are available in 
reference 3 and indicate that the potential core extends 
approximately two jet diameters from the nozzle exit for all jet 
velocities. As a test of the sensitivity of predicted surface 
pressures to potential core length, a parametric study using 
potential core lengths of 1.5 and 3.0 jet diameters was carried 
out for Vj/V w » 3.43. Results of this study indicate nearly 
identical predicted pressures using both potential core lengths. 
Since flat plate data have shown (ref. 7) that a jet in a 
crossflow has a shorter potential core length than a jet under 
static conditions, 1.5 jet diameters was selected as the 
potential core length for all calculations included in this 
report . 

A series of expansion rate curves corresponding to a 
potential core length of 1.5 jet diameters is shown in figure 
21. The Vj/V^ = 1.96 and 3.43 curves were obtained by 
extrapolating the Vj/V w =4, 6, and 8 expansion curves 
using information from figure 22. This figure presents (t/D) 
vs. Vj/V^ where (t/D) max is the value of t/D at S a /R = 20. 

This curve is utilized in the same fashion as outlined in 
reference 18. 


Correlation .- The purpose of the correlation for the body of 
revolution case is the same as that for the flat plate case; 
that is, to isolate the viscous or non-potential effects on the 
jet on the body. It is expected that the predicted potential 
pressure distribution away from and upstream of the jet will be 
in good agreement with experiment. The predicted pressures near 
the jet and downstream of the jet will likely be in poor 
agreement with experiment due to viscous effects. Equations (6) 
and (7), which apply to the flat plate correlation, are also 
applicable to the body of revolution correlation. 
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In a manner similar to that done for the flat plate case, 
(AC ^viscous va * ues are obtained by comparing theoretical 


pressure coefficients with experimental data at chosen $ and 
x/D stations for each Vj/V w . A data base is set up which 
consists of a (Ac^) v ^ scous array for each jet velocity ratio. 
Linear interpolation is used to determine the correlation 
factor at any Vj/V M , §, and x/D values for which correlation 
curves have not been determined. 


RESULTS 


Flat Plate 

t 

Correlation curves are presented in figures 23, 24, 25, and 
26 for jet inclination angles 15°, 30°, 4S P , and 60°, 
respectively. The number of jet velocity ratios for 
which correlation curves are presented for 0 = 45° and 60° 
are limited by the available data. The correlation values were 
determined for most velocity ratios at $ = 0°, 30°, 60°, 90°, 
120°, 150°, and 180° and r/D = .75, 1.0, 1.5, 2.0, 3.0, and 5.0. 
Correlation curves obtained using data from reference 6 only do 
not contain correlation factors for r/D = .75 and 1.0 due to the 
absence of data at these radial positions. Also, correlation 
factors for all r/D .values could not be obtained at each 0 value 
for the data from reference 7 due to lack of data cit various r/D 
positions. For cases in which two independent sets of data were 
used to obtain correlation curves, the given velocity ratio is 
the average of the velocity ratios from the two set/? of data and 
each correlation value is the average of the correlation values 
obtained using each data set. In most cases, multiple data sets 
were not available and the correlation curves were obtained 
using only one independent set of data (refs. 6, 7, or 28). 
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The results shown in figures 23 through 26 are similar in 
form to those presented in reference 18 for jets exhausting 
normal to a flat plate. The (Ac p ) viscous values are usually 
smallest in the region upstream and away from the jet and 
largest in the region immediately downstream of the jet and in 
close proximity to the jet, as would be expected. It is noted 
that the erratic behavior of the r/D « .75 correlation curves 
for the region 90° £ 6 1 150° is very similar to that 
seen in the 0 ~ 0° case in reference 18. As was previously 
discussed, the correlating factors are maintained in a data base 
form as shown in figures 23 through 26. Any required factor can 
be obtained by interpolating between curves for a given 
velocity ratio and inclination angle combination. 

The pressure prediction method, made up of the blockage 
model, the entrainment model, and the viscous correction 
factors, has been applied to jet velocity ratio and inclination 
angle cases for which two independent sets of data are 
available. Since these data are the same data used to obtain 
the correlation factors, the purpose of these results is not 
strictly to show comparisons of predicted and experimental 
results, but more so to show how close the predicted values are 
to the data relative to differences between the two sets of 
data. Obviously, it would be more useful to compare the 
prediction method with independent data, but adequate additional 
data are not available at the present time. 

Data comparisons of Cp vs. r/D at 0 = constant are presented 
in figures 27 and 28 for the following jet velocity ratio and 
inclination angle combinations: 3.9/15°, 5.0/15°, 4.0/30°, 

and 8.15/30°. Data are not available at r/D = .75 and 1.0 
from reference 6, and therefore the correlation factors at 
these radial positions reflect the data of reference 28 only. 
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As discussed in reference 18, the success of the prediction 
method is dependent on how well the various sets of data agree 
with one another. In regions where the data set are in good 
agreement with one another, such as V./V =3.9, 0 *15°, 0 * 60° 

J 00 

[fig. 27(a)], the predicted results agree well with both sets 
of experimental results. In regions where the data sets are not 
in good agreement, such as V^/V^ - 5.0, 0 » 15°, 0 « 30° [fig. 
27(b)], the predicted results show an average of the two sets 
of data. Perhaps the most significant aspect of these results 
1» i 3 that for 0 > 90°, the viscous and inviscid contributions are 
of comparable size (see, for example, figs. 24(b) and 28(a)) 
and the separation of the phenomena into a predicted potential 
and an empirical viscous correction results in uniformly good 
results over a wide parameter range for engineering calcula- 
tion purposes . 

Pressure distribution data on a flat plate with a jet 
exhausting normal to the plate were also obtained in the 
NASA/ ARC tests of reference 28. These data, in conjunction with 
data previously used to obtain the correlation factors presented 
in reference 18, were used to improve the viscous correction 
factors at 0 =0°. These correlation curves are presented 
in figure 29 for jet velocity ratios 2.2, 3.9, 5.0, 6.1, 7.0, 

8.0, and 10.0. New correlation curves for Vj/V = 1.0, 

J 00 

1.67, and 3.33 were not obtained, since the new experimental 
results did not include these velocity ratios. 

Comparisons of the NASA/ ARC data and predicted results 
obtained with and without the NASA/ARC data included in the 
viscous correction factors are shown in figure 30. Comparison 
of the new correlation curves with those presented in reference 
18 indicates that the addition of the NASA/ ARC data, in general, 
had little effect on the correlation curves for 0=0°. The 
largest effects are Seen for r/D _< 1.0, where experimental 
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pressure data are difficult to obtain due to large unsteady 
effects near the jet. It is noted that correlation curves had 
not been previously obtained for Vj/V^ = 7.0; therefore, 
predicted results presented in figure 30(e) are available only 
with the NASA/ ARC data included. 


Body of Revolution 

To verify that the potential model results give good 
agreement with data in the region of the body of revolution 
where viscous effects are small, data comparisons of surface 
pressures were made for Vj/V M = 1.96 and 3.43. These 
velocity ratios are the- lowest and highest velocity ratios for 
which data are presently available (ref. 3). 

Comparisons of measured and predicted surface pressures for 
Vj/V^ = 1.96 are shown in figures 31 and 32. The predicted 
results at this velocity ratio were obtained using the blockage 
model alone, since Yeti' s entrainment model is singular when 
Vj/V^ < 2.35 (see refs. 14 and 18 for details). Figure 31 
presents axial pressure distributions for several constant 
positions, and figure 32 presents circumferential pressure 
distributions at several axial stations. These comparisons 
indicate good agreement between theory and data in the regions 
upstream of the jet ( x/D >_ 0), although the jump in experimental 
pressures between x/D = 7.0 and 6.0 is not shown by the 
predicted pressures. This may be due to separation at the 
cone/cylinder junction ( x/D = 12.9), which cannot be predicted 
using the NSRDC program. In the strong viscous region 
downstream of the jet, predicted pressures compare poorly with 
measured pressures, as would be expected. In regions aft of the 
jet where viscous effects, and jet effects in general, are small 
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(x/D < -5), measured and predicted pressures are in good 
agreement. 

Predicted pressures for V./V M = 3.43, obtained using both 

J 

the blockage and entrainment models, are compared with measured 
pressures in figures 33 and 34. These comparisons indicate the 
same general trends seen for Vj/V M = 1.96, although 
comparisons between predicted and measured pressures at x/D = 

1.0 and 1.5 are not as good as those seen for Vj/V^ = 1.96. 

This may be due to inaccuracies associated with the use of a 
flat plate entrainment model to predict entrainment effects on a 
body of revolution. For the flat plate case, it was found that 
entrainment effects on surface pressures become important at 
large jet velocity ratios. This is most likely true for the 
body of revolution case as well. Lack of body of revolution 
data at higher velocity ratios prevents further investigation of 
the adequacy of the flat plate entrainment model for use with a 
body of revolution. 

Correlation curves are presented in figures 35 and 36 for 

jet velocity ratios 1.96 and 3.43, respectively. The 

correlation values were determined at a range of <p values from 

0° to 25° and a range of x/D values from +1.0 to +7.0, 

where "+" and represent positions upstream and downstream of 

the jet, respectively. Since predicted pressures could not be 

obtained at tfe = 0°, (AC), . values for this angle were 

p viscous 

obtained by extrapolating the theoretical results to <j> = 0° 
at each x/D position. For I x/D I > 7 and <j> > 25°, ( AC p ) v ^ scous 
has been set equal to 0.0, since jet effects are negligible in 
these regions. Also, since data are not available in the region 
I x/D | < 1, ( ACp ) v ^ scous values are obtained by interpolation 
for x/D values in this region. 
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values are 


As these correlation curves show, (Ac 

p VXoCUUb 

small in the regions upstream and away from the jet and large in 
the strong viscous region immediately downstream of the jet. 
Also, the large correlation factors decrease with increasing <j>, 
as would be expected. 

Since these correlation factors have been obtained using 
data from only one source (ref. 3), these factors are not used 
here to calculate new pressure distributions on the cone- 
cylinder model. It would be useful to compare predicted results 
with independent data, but adequate additional data are not 
available at the present time. 


CONCLUSIONS 

A correlation method to predict surface pressures on a flat 
plate in a crossflow with a jet exhausting normal to the surface 
was extended to jets exhausting at an angle to a plate and to 
jets exhausting normal to the surface of a body of revolution. 

As in the previous method, the correlation values of 
viscous-induced pressure coefficients are defined as the 
difference between pressures predicted by a potential model and 
those obtained from experimental data. 

For the flat plate case, blockage and entrainment effects 
are accounted for using modified versions of the vortex 
quadrilateral panel model and the sink distribution model, 
respectively, of reference 18. The viscous effects of the jet 
are represented by the correlation values. 

Comparisons of measured and predicted plate pressures in the 
vicinity of inclined jets over a wide range of jet velocity 
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ratios and inclination angles are generally good. Correlation 
values for certain ranges of jet velocity ratio and inclination 
angle are to be considered preliminary due to inconsistencies in 
the available data (ref. 28) for those ranges. 

For the body of revolution case, mutual interference effects 
between the body and the jet surface are modeled using a source 
panel method. Entrainment effects are accounted for using the 
flat plate sink distribution model due to lack of an available 
method for a body of revolution. The viscous effects are once 
again represented by the correlation values. 

Predicted surface pressures, which were obtained for 
Vj/V^ = 1.96 and 3.43 using the potential method, were 
compared with measured pressures in the vicinity of the jet and 
showed good agreement for regions where viscous effects are 
small. Lack of available data made it impossible to assess the 
method for jet velocity ratios greater than 4 and jet-to-body 
diameter ratios greater than 0.1. 

In addition to extending the normal- jet correlation method 
in the aforementioned cases, new normal- jet data were compared 
with results obtained using the normal- jet empirical prediction 
method with and without the new data included in the viscous 
correlation factors. In general, the new data had little effect 
on the predicted results, which lends confidence to the use of 
the method for normal- jet cases. 

In summary, the following comments are made regarding the 
applicability and use of the viscous correlation prediction 
method. First, the correlation factors for subsonic jets 
exhausting normal to an infinite flat plate can be used with 
confidence. These factors were obtained using data over a wide 
range of flow parameters and from several sources. Data 
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comparisons indicate that differences between predicted and 
measured pressures are usually less than the differences between 
different sets of data. This indicates that the correlation 
method should give accurate flat plate pressure distributions 
for the normal- jet case. For jets exhausting at an angle to a 
flat plate, the correlation factors for low and medium jet 
velocity ratios (V^/V _< 7) and jet inclination angles (0) 

less than 45° can be used with confidence. Inconsistencies 
in the data for high velocity ratio jets suggest that the 
prediction method be used with caution for Vj/\ T ro >_ 7.0. For 
jets exhausting from the surface of a body of revolution, the 
correlation method is to be considered strictly preliminary. 

This is due to the lack of available data, which limits the 
applicability of the present method to jet velocity ratios 
less than 4.0 and jet-to-body diameter ratios less than 0.1. 
Additional data are needed to extend the method to higher jet 
velocity ratios and to larger jet-to-body diameter ratios 
typical of full-scale fan-in-fuselage configurations. 


RECOMMENDAT IONS 

It is clear that the prediction method described in this 
report (and in reference 18) depends heavily on experimental 
data for determining the viscous correlations, evaluating the 
potential model, and assessing the accuracy of the overall 
prediction method. In many cases, the available data are very 
limited, and the recommendations have to do principally with the 
need for further data. 

For the flat plate case, it is encouraging that inclusion of 
the most recent NASA data for normal jets has little effect on 
the viscous correlation developed previously. The principal 
need at this time is for additional data for jets inclined at 
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angles other than 90 degrees to the plate. Two types of data 
are required. First, data on the jet properties (centerline 
location, expansion rate and entrainment) over a range of jet 
inclination angles and velocity ratios are required to assess 
the accuracy of the jet model and its ability to predict surface 
pressures upstream of the jet. Secondly, data on surface 
pressures are required both to add to the data base for deducing 
viscous correlations and to independently check the overall 
method . 

The body of revolution case has been investigated much less 
extensively, and the two types of data mentioned above for the 
flat plate are required also for bodies. The body case also 
introduces more geometry parameters, in the form of hole-to-body 
diameter ratio and shape of the body (cross-section and axial 
changes) . The most important need is data on configurations 
with larger diameter ratios than are currently available. 

Finally, as noted in the previous report (ref. 18), 
measurement of jet turbulence level and mean velocity profile, 
in addition to induced surface pressure measurements, would be 
desirable. Such data could aid in the understanding of real jet 
effects and in the improvement of the present prediction method. 
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Figure 6.- Measured potential core length for a jet exhausting 
from a flat plate into a subsonic crossflow. 
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(a) Vj/V^ “ 4.0 

Figure 11.- Comparison of measured and predicted 
pressures on a flat plate from which a jet 
is exhausting into a crossflow, 6 =» 45°. 
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Figure 11.- Continued. 
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Figure 13 . - Comparison of measured and predicted 
pressures on a flat plate from which a jet is 
I issuing, V./V = 6.67, 6 = 15°, p = 60°. 
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(a) Reference 4 data. 


Figure 15 Effect of jet velocity ratio on measured 
pressure distribution on a flat plate from which 
a jet is exhausting into a crossflow, (3 = 0°. 
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(b) Reference 6 data. 
Figure 15.- Continued. 
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(o) Reference 28 data, continueti 
Figure 15.- Continued. 
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Figure 22.- Maximum jet thickness variation with jet 
velocity ratio for a jet exhausting normal to 
the surface of a body of revolution, S /R 13 20. 
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Figure 23.- Continued. 
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Figure 24. - Continued 
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Figure 24 .- Concluded. 
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Figure 25 Correlation factors for viscous portion 
of the pressure coefficient induced on a flat plate 
by a jet exhausting into a crossflow, 0 = 45°. 
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Figure 25.- Concluded 
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Ficmre 26.- Correlation factors for viscous portion 
of the pressure coefficient induced on a flat plate 
bv a jet exhausting into a crossflow, 

Y J v./ v = 12.0, 0 = 60°. 
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(a) Vj/Vo = 3.9, continued. 
Figure 27.- Continued. 
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Figure 28. ~ Comparison of mecisured and predicted 
pressures on a flat plate from which a jet is 
exhausting into a crossflow* 8 = 30°. 
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Figure 28.- Continued. 
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Figure 28.- Continued. 
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Figure 29,-Correlation factor for viscous portion of the 
pressure-coefficient induced on a flat plate by a jet 
exhausting into a crossflow, 6 - 0°. 
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Figure 29.- Continued. 
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Figure 29. 









Figure 29.- Concluded 
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Figure 30.- Continued. 
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Figure 30.- Continued. 
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Figure 31. 








Figure 31.- Continued 
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Figure 32.- Continued 
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Figure 34. - Concluded. 



Figure 35.- Correlation factors for viscous portion of the 
pressure coefficient induced on a body of revolution bv 
a jet exhausting into a crossflow, Vj/V = 1.96, 0 = 0°. 
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Figure 35.- Concluded. 
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Figure 36.- Concluded. 





